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ABSTRACT 


For certain assumptions, an analysis of multiples HAD 
signals results іл a reasonable estimates for the localiza- 
Rot а target relative <c thes MAD platform. This is 
achieved by using selective approximations to linsarizs an 
initially nonlinear problen. The simulation igneres noise 
and requires an =stimate of the nagnitude of the target 
magnetic moment components. Results indicate shat she be 
localization estimates are achieved when the platform 
cardinal headings, and when the <target moment has a strong 


vertical ccmponent. 
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I. INTRODUCTION 


А. BACKGROUND 


Pes = псе "the operational inzfoduction of the MAD 


(Magnetic Anomaly Detecticn) зүхэл in 1944, i+ has реет а 
single sensor system limited “o iadicating  de-sctior only. 
The short ranges of the early varsions allowed a general 
localization if a signal was detected. The use or nultivie 
sensors should allow for a good 3st салса 


imate of the 
location relative to the sensor platform, assuming that the 
arrangement of che sensors is such that the received sigqnals 


are linearly independent cf one azother. Previous attempts 


h 


h 
to solv2 this problem have in general focused on the use of 
only two sensors, but did not develop a usable solution t 
Seecatsonal completion. Work has been dene on this topic by 
Wynn et al., [Ref. 1], using г superconducting aradiometer 
array. The measured signais in thair work were subjected to 


a novel signal processing techni чаше и за “Re. Labor 


py 
| 


tory, was used to invert the dipole fieid equation to deter- 
nine both the position and the monent vectors of a dipole 
signal source. This thesis approaches the same problem by 
simulatirg five standard total field maasurements. As shown 
20 Figure 1.1, four magnetometers moving in parallel paths 
are used to obtain these neasurements. The advantage of this 
Roach compared Зо that of Wynn 2t al. is that it utilizes 
E coro *techpolsgy in conjunction with approximations, 
delineated in Chapter 3, БОК == 29 produce a usable 
suilt. 





B. APPROACH 


ШПЕЕ Шепа Е localizing а -arget from any platicza 
myoives at isasc three unknowns: Су» = алам 
ments Of he position vector from the platform то the 


target. Іп the case of using MAD, three additional unknowns 
ғ 


(t 
Er 
(р 
«+ 
= 
2 
D 
сі 
(Л 
ЕЗ 
NY 
ч 
m 
р 
(а 


m and в the components 


х”? Nye т, 
moment, are also involved. The Tas: 


= о 


unknowns that must be solved for in 
Serge: using MAD. білсе з normal 512 
be approximated by the magnetic field 
iż was decided to use thea dipole 2а 

fRef. 3] +o generate the signals us 
БА 00 1.1 3escribes :he MAD signal, and w-11 
mou the dipole $quation in chapter 2. Equatlo 


Me (tet matte) 76743 (act meee) (re+sa+re)}/>S (eqn 1.1) 


Цеа 17) the unknowns. The posicion vector ina equation 1.1 


Ae rined as aqong arom the target zo the platticrm: 5115 is 
Меке =-аџегс= cF using zhe ёїроі= equation, and will bs 
Ер ал зад later. 

Ths cartesian ccordinates аге 32:ined I n Figuss 1.2. 
О Ол ест тор of platicim motion, anc “да y 
ENDS -5 үегс1сахл 30 apen Ort о т д 


IE тестот, зи 2,0 and 5, , 246 assunsd co he know 
ЖЕ сз се 2: 15 very iircz2cult <3 explicitly solve a set 
or ron-linea 


Meeta] Casas., The f25st case 1S 2 1211 
Epsroet!s position is assumed known. For this case, =quaticn 
fet 15 line Deco tires саклама п, п, псп If the 


E 
гэе measurenents are linearly inasvoandent, che 


RulS, [ Ref. 47], may be used to 50173 for these unknowns. The 


10 





second case is the mere difficult one, and it is the orcpien 
investigated in this thesis. In this case, the assumption is 
made that the target's magnetic moment is known. = Rake 
E-ducss equation 1.1 Into a non-linear equation in only 
three unknowns, г,, ry, аллаа = ле клада Version 
песама -2оп 1.1 could be manipulat= into a linear form by 
using several sensors and appropriate ccmbinations of 
Signals, the positicn of the target could sesily be solved 
БОТ. 

This thesis will develop an approximate solution то this 
Beoblem by using combinations of he total 21814 sicnels of 


an idealized multiple MAD systen. 


11 
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Pigure 1.1 


20125. 


Spatial Relationship of Sensor Posit: 
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13521 
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45. 


2434 
$9 
мере р 

144 40, 3714 


153 





II. THEORY 
The total field magnetometer measures the aagniczude of 
the vector sum of the ambient magnetic field of the earth 
and the magnetic field of а magnetic dipole moment. The 
peta of the dipole will used то represent the magnetic 
field of the target. The magnitud= of the vector sum is 


called the total field as given by 


5» ES == 
8, = [Be + byl 

сересеп the earch’ magnetic field vec-or and 54 
represents the field dus to a magnetic dipole. бы со 
generali case ol interest here, the magnitude of the dipole 
field is much less than -h2 magnitude of the sarth's field. 
Under certain conditions, а useful appreximation со “із 
total field may be derived from the geometry shown in Figure 
201 . 





In =his figur2 a very small vector is added to a very large 
> т? : 

опа. Вы о Е Че 2.1 аса пот то scale. by -5 60 

be considered several orders of magnitude smaller than Be 

The application of basic geometric formulas for the sides of 


a triangle and tne angle between sides yields 


Каже - 5-2 € ha? - 2bB,cos 


14 





ош тіле сеопе=су о: Е:ди== 2.1 - Гал terms cf e this squarion 


becones 


BL? 7 Bg? * bj? + 2b4B,cose 


The B,? term can be factored out of che right hand пора = ос 


the equation to produce 
В, 2 = Ве? {1 + (54/8,)2 « 2(Ь4/8.)с>59) 


Taking the square root yiélds equation 2.1. 
B, =B, {[ 1+ (by/Be) 2+2 (bA/Be) cose ]¥7} (equ 291) 


Since by assumption by is much 1355 than Be, “ке term 
EE is mich 1555 than  2(54/B.)cos8 and can be consid- 


eced negligiblas. This reduces equation 2.1 to 
8, =8, ([ 1*2 (b4/ B.) cose] v?) (247222) 


The general expansion formula given by 
n . 
Ш) = ТЕ nx * n(n-1) x 2/2! * Higher Order Гегт< 


БІ Пе арр1161 ко залгтісп 2.1. Poraba са 
ВОСИЕВ „Баије +2 (54/8;) 2059. Substituting chess values into 


the general expansion formula gives =quation 2.1 as 


ШО = Вот “ .5{ (by /Be) 2 + 2(bg/Be) С259 


г 172 (5,054) t (254/78,) cos8 J^ * H.-O.T.J 


t 


eeng the appropriate <2rm and multiplying through by the 


ee ficients results ir 


БЕРУ Гы.) =] + (ра/8е) 2058 - 1/8[ 5478. ]* 


- 4/8 b, /B, ]3cos9 - 4/8[ ,/3, |2со52Ө + H.0.T.} 


Miemeub=d and higher order terms will be considered +o be 
negligible. Dropping 21952 terns апо из = алата th 


remaining ones yields the forn 


15 





БИН ИООФ соген — -5(54/В-) 2с0528 + .5(54/56) 2) 


t 
Using the relation 511029 + со5згд = 1, this raduces te 
8, Ве (1 + (Ба/Вь ) с059 + .5 (by /Be) 251220) 


Миш ер јула through by Bo, the second term can be 
— 
Ep tren in terms of the dot product ba * Be /Be - Пи = 15 
-— A : Е : = 
ас осеб со byes if the unit vector $ is defined as Be / Pe, 


and the resulting equation is 


8, - B. + Dye * .5(b4/8e ) ?sin?e 


се ба 15 on the order 5f InT while B& is on the orde- of 
ООО Вт, the error tsrm involving the sine is nagiigible. 
For the stated assumptions, equation 2.1 can be rzeoresented 
5у 
т A 
8, и ЛЫС (297762551) 

nnus thesis, sd, will be assumed constant. In addiction, 
ams simulation will assume that ro noise exists. Thus, the 


tal field magnetometer will 21635113 the anomaly in the 


earth's field defined by 


= 


S(T) = he 


Т» 


In this expression, by represents žhe m 


Q 
= 
(D 
(4 
ps 
Q 
Fh 
- 
(Г, 
1-4 
ku 
О 
tt 
fu 


а 
Боје. by -5 defined [Ref. 3) by the =quazicn 
== - 
ШІ (2) - -п/с3 * 3(Wer)r/rs (eqn 2.4) 
m is the magnetic dirole moment in units of aTft3. F is the 
position vector to the field point and is measured in feet. 
о 


НО is much less than Be, then =quation 2.3 ала ғала-зі 


2.4 result in the equation 


16 





S(F)=- (me2)/r3+3 (Mel) (2»2)/с5 


(сап 2.5) 


Miss equation 1s an alternative basis for the derivation of 
Ене Anderson's Functions, found іп [Ref. 5]. The Anderson 
Functions area commonly used approximation for the signal 
received by a MAD sensor. ло ccoupmeer simulation in this 
eesis is passi on equation 2.5 directly, since it is ina 

tation which is easier *c uss than ће Arderson's 


PED ESen 
n 


representation of a MAD signal. 
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Bhucethesrs Wall develop and demonstirate 22 approx 
Бион chat generaras a гвазолавіев solution cf tine target's 
position relative to a platt 
"measure" “he magnetic signal as the — 
traigh iine encounter nea 
аз follows: 

ОЕ а: роти equation is a valid T[epresentalicn of 


mre Magnetic field of the target. 


as 
a 
153 


E Dy is several orders of magnitude ca 


a + 
Sm Ucing only s -raigh line en- 


(0 
Си 


.) Measurements ar 
anters. 
= -==<0е-'<= magnetic dinsie moment is a2 Хдомг 
тесто 


КОШТ е Басын magnetic field vector is Know 
2 


ута НА 


Ма тебе об агу ksnd (opv--onmss*t, сепзсг, 


“) 


т 

Пт п of the earths 21214 15 
пе 

с) 


2.) Ihe measurements are criented along four lings 
эс го DDLen zed wich respect to cae flight path as 


shown in Figurs 1.1. 


Бле ths dipola equation is written with tha moment at “hs 

йіп ofa cocrdinats system, tus vcalculsctOoNS oS таз 

Sana fislid values usec in the simulation are made wich the 

mame: locatz3 a- the origin. AS Will be seen, he progran 

BPepecrms 2 COOrdinate transformation in crder to produce the 

Meet Of zarjet posicion relative “o the platform. АЁ^әг 
O 5 


Ш COOrdinate transformati 


meso: Platform and meves wich iż. 


ct 
(D 


The five nsasurements are maie along four separa: 
В 


Ф 
|: 
(3 
«4 
о 
|--3 
«ч 
(D 
{u 
• 


Шап: paths, so only four di££ezrzsn: "sensors" a 


13 





Sensor 4 could бе герјасга by 22012124 а time delay in the 
measurement of sensor 1. meS nals Generated аз thess 
measurement points are measured at sach time step. The posi- 


tion vectors can all be described relative to the origin 
ЕО d= che target, in terms of the position vector c 
sensor 1. The measured signals ar2 represanted by the 


following equations: 


5(2\)=- { (щазтая це) (г1)2-3 ( цх+цу+ ЦЕ} (хэ4 уз 23.) } /г 15 

Меј == (пеел= пе (12) 2-3 (дх+ц[у-а ] +32) (к э+Ёу-9 Jerzq) } /= 25 
5 (73)=- { (metas mg) (3) 2-3 (ај х-а ]*ay + д2) (L Хх-216%у<%25))/235 
S(5) -- ( (matus ag) (r0) 2-3( xe nyen[ z-d]) (xy sp( 273 ]9) ) /: u5 
S (55) =- { (matne+ mg) (т5)^=3 (ај хча Jay +mz) (( хе јећђус 25) } /= 55 


ееисогтегпјепса, the following dummy variables are used in 
the computar program anrd are dsfinesl as: 


G = П,2х%1/3,%т,2, 

V1 - 3ш,6, - G 

V2 = 3 (m =, + mye ) 
V3 = За,“ - G 

74 = 3 (n, 2 + R2) 
ШИ - 31-2. - С 

М 3(ш е, - m 
51 = s (3%) 

БИ - S (52) 

53 = $(=3) 

54 = 5(=5ђ) 

55 = s(d) 


Tha use of “these variables reduces the above equations to 
Пе с сы па form, which are the j2tining equations for the 


mals in he simulation. 


S1= (V1x2+V2xy+V3y2+ Yüyz+V5 z2+V6xz) /r 15 
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S22 (V1x?*V2xy-V2xd*V3y2-2V 3yd*V3d? «TV 4у2- 7442 +752 2 + 9 6Х 2) /525 
S3= (11x2-2V1xd+V1d2+V2xy-VW2dy+V3y2+V4yz+V5zZ2+V6xz-V612)/7353 
6Цт-(У1х2%У2ху%У3у2%Ү4у2-У5у4%7/522-2У524%Ү542 чуб ха -Убха4) /-45 
55- (У1х2%2У1Х24%У1424У2ХуУ%У24у%У3ү2%У4Уу2%У5224У6х24%Ү642)/255 


The unknuncuns of interest are x, y, and z, which are related 
О 


mn hbe рог1:1оп 7225 


Е1 = х% + уј + zk (ean 3.1) 
ro = xt + (1-1) 31 + Zk (салт 222) 
23 = (х-4)2 + yf * zk (ғап 3.3) 
Fh = xt + yf + (z-djk (eqn 3.4) 
25 = (х+а)т + 77 + 2К (can 3:5) 


Шо нсаре He =preplem it will be assumed that the 
 ШБЕпапатсгс 22 һа equations for 51, 52, 53, 54, and S5 
тег above ere all squel to ст, 1.3., 

Ш = 11 = 12 = 13 = 54 = 25 
Пи - апа-:ов5 гог 51, 52, 53, 54, and S5 wirh +: 
Л Э 


раното Со. Со 


$2-$3= { (12-291) х+ (243-12) 7+ (14-16) 2+ (91-\3) 4} (-а/=5) 
52-55={(12+271) х + (213 +12) у+ (Үч+ү5) 2+ (1-3) 1) (-4/г5) 
53-51-12У1х%У2у%У652-714) (-4/-5) 

ЕЕЕ 5 |І-(-2У/1х-У2У-762-УЛ1а) (-4/-5) 

б -- <= ((гУ1– Уб) х + (У2—7%) уч (96-2У5) 74 (Т5-У1) 41 (-d/rs) 
355-S42 (- (2U1* T6) x - (V2«VU) y- (V6*275) ze (V5-V 1) 3) (-4/75) 


Since +) 


(D 


diusmmnafFPo-c are assumed co hs i¡centical, if these 


six equations are divided in pairs о сп 
вета tne non-linear гә егт ani the quantity d are 
canceled ояз. ВЫ ЕЕ в ты» 201104: 84 sauations chat 


define cha пем variables Al, A2, ani A3 as 


= (52-53)/(52-55) (eqn 3.9) 
22 = (53-51)/ (55-31) (eqn 3.7) 
АЗ = (53-54) / (55-54) (247 3.8) 
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А се ом су усзгізбізе for convenisncs, 


C1 = A1(V2+2V1) - (V2-2V1) 
C2 = 11(273+72) - (2У3-72) 
C3 с А1ЦУ4%У6) - (УЦ-76) 
Шри (Џек - 11(/1-У3) 1а 
СЦ = А2 (241) + 241 

C5 = A2(V2) + V2 

C6 - A2(V6) * V6 

K2 [V1 - A2(V1) ]d 

C7 2 A3(V6*2V1) - (V6-2V1) 
С8 з АЗ(Ү4472) - (14-42) 

C9 = 43(2У5%У6) - (2V5-V6) 
КЗ = ((71-75) - A3(V1-V5) Ja 


ENDE cos 3.5 - 3.8 can be wr2t-er in che following form as 


ex + C2y + C32 * Kl (eqn 3.9) 
Sox + C5y + Coz = K2 (eqn 3. 10) 
mix ж CSy + C9Z = КЗ (саг 3.11) 


Miese тіссе Equations are linear in 52245 of х, у, and z. IN 
Ба слог, they are linearly inde 
the program in Appendix B. А sample c 


linear independence is shown in Ap 


гол 0 
Р 
о 
fv 


linearly independent, кешек "е Ко 


1 
x 
2 
Кола 3.9 = 3.11 ог X, у, 2112. The co 
in Appendix A was writtsr in n 
program assumes that the magnitude 22: 5: 
Moment is 5x198 ПТ“ЕЗ. Are 
buus2ng the vaiues listed by Fromm in [Ref. 2], who iists a 
Ament of 108 to 2x109 cgs units for a submarine. One cgs 
ИЕ 115 appro ximateiy 3.35 11253. Th 
program is 2 rounded out average of Fromm's values. 
Computations are done using che s 


meme. di: pole at tae origin. However, as indicated sarliec, 
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PRE NOULDU: 2s corrected by a coordinate 
ас 11 placform is ar the oràigir, 


Ncogriven with respect to the platform. 
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о 000 2000 3000 4000 5000 8060 7000 8600 9000 10000 
DISTANCE ALONG TRACK ———> 


LEGEND 


o = CALCULATED 
à = ACTUAL 


Figure 4.2 Initial Run, y position. 

Prom the orisntation of the coordina ш shown in 

Bene 1.2 , 2 1S in the direction of flight, у is straight 
о 


ШИЕ) апа х 55 off the tight side of the olatform. Since 
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Z POSITION 
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21111 
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о = CALCULATED 
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Figure 4.3 їп1*1а1 Воп, z position. 
this thesis investigates straight line anc 


С 
ос а!) and y(actual) cocrdinates will тэта? 
n 


EucuUughout any run for thes coordinate oris 
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Figure 4.4 


CORRECTED Y POSITION 


2000 3000 4000 5000 0800 7000 3000 3000 10000 
DISTANCE ALONG TRACK -——-> 


LEGEND 


о = CORR. CALCULATED 
s = ACTUAL 


Initial Run, corzec:ed y position. 


11 therefoza be nađe wiza their respec- 
ТОО eee” aS сопсталсс. Таз 2 coordinato will 
ШЕП ат ОЕ Ene 201:141 Dass j=sczibed, 
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actual) would be constant et zero, y(actual) would De 
constant at 1000 feet, while z(actual) would pregTsss 2758 
КОО feat 075000” Безт. Figures 4,1 - 4,3 show 2 compar- 
Gson of the actual and calculated targst positions as the 
PeseLOrm cOmplet ss the initial Fun described earl 

will be noted that the calculated and actual x, у, аа? 
puc ons are almost identical at the start of the run. This 
femme coincicence. "Crossover" points such as this occur 20: 
all їїсєв ccordinates simultansously during a rus, failing 
Шогег ог fureher away from CPA, depending upon the choice 


of parameters. 


Un chis inicial TUL, Ме ле Вов“ о the "worst" 
Bee tates occur tor the y coordinata. This is less signifi- 
cant than it seems at first glancs, however, ас the y cocr- 
dinate is the vertical separation between the platform and 
НӨ target. If the target is submerged, this value is 
Питом, DUT tne altituds of the platficrm is known. The 
те can be no closer to ithe piatform than this altitude, 
Since he target is physically unable to rise above the 


V 
ocean surface, As an example, if a tarzgs- Gspth of 200 “зат 
Ї C 


is assumsd, the initial y separation of 1000 


ШИРЕ псе Zn this Simulation to be 390 f¢2t. Using this чаіца 
Шаға bourd oa y, it can be ss¢n in цэг у 


calculat=d position is somewhat impro 
ms that the particular arrangensnats 
штите 1.1 des show promise in n 


E cecmetry of Multiple MAD “for acc 


purposes. Меке "лошу, чатшасіспаз of? paranmnekszs were +2224 
от дет *o de-ermine wha- effect changing thea would have. 
in all cases, parameters not specifically mentioned іг а 
Eon Well retain the values as for the inicial run 
described above. Деш па о = а: у ке behavior occuring at 
дет саз Places on the graphs are a result of the deterninant 
ENECOG denominator in the Cramer's Rule subroutine having a 


value close то 2350. 
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В. HEADING CHANGES 


сэ. the localization process will be inves- 
шдатед аза Euncticn of heading. Compared to the initial 
conditions described above, the crossover points w 
to be moved further cut past 5009 fəs: prior t 
heading of 09. While a specific example of + > not 
shown, the behavior cf th2 crossover points for the platform 
headings of 309 and 60° may be observed by comparing Figures 
4.1 - 4.4 to the respective portions f Pigura 4.5 . The x 
position estimates were better for а North-South heading 
than for any other platform heading, although the estimates 
on East-West headings were also very good. For a 90° heading 
there are two crossover points located equidistanz at 1200 
feet from CPA. Calculated positions on all headinas ara 
generally gcod, excepting those around 600, where the x 
estimates snows a larg2 discrepancy at distances greater 
than 3000 feet prior to CPA. The y 2stimates are nor Ss 
icantly affected, but the z estimates are good only within + 
1200 feet of CPA. Nevertheless, sven in this worst-case 
Situation,  r-z2sasonable calcula-:ed coordinates are available 
оса 1221101 2iong a Significant pozticn of the track. 

As shown in Figures 4,5 - 1.8, headings in the quadrant 
09 - 909 reasonably describe what happens on any heading, 3s 
the graphs for analogeus headings (509, 1209, 2409, and 2000 
in this example) have similar shapes, although «he crienita- 
tion may differ between quadrants. Diese sinilarities arə 
apparently dus tc  *he fact tha* analogous heading in quad- 
pars | and 3 and quedrants 2 ас24 4 за in fac 
headings, and so one would expsct -hs graphs + 
same shape; however, the ending point of 


ШЕ стагтіпа ролпт of the other. Iha re 


ЕЭ 
f» 
(1 
>: 
O 
"2 
Л 
=. 
! 34 

о "а 
СУ 
iD 
ГЕ 


adjacent guadrants probably exists for similar reas 
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С. TARGET DIPOLE MOMENT CHANGES 


ro ehis Section, she horizontal component of was 
oriented in the 00 (magnetic North) direction vice 355%, and 
the vertical component of ПШ was varied. The best positioning 
was obtained for RH oriented vertically up,i.é., no hori- 
zontal componen®, and an example for that case is shown in 
figure 4.9 . Poor positioning was obtained for tne vertical 
component of 3 pointing between thea horizontal and 45% up. 
Mne worst case shown is Figure 4.10, whers 2 is oriented 
мо zoteally pointing dua north. With the notent oriented at 
an angle below the horizontal, reasonable positioning esti- 
mates were cbtained. There was less variability in posi- 
EUonong for moments pointing below the horizoztai than for 
Enose po2n*ing above the horizontal. 

The results for varzations ia the orientation of the 
target dipole moment suggests that tha positioning 
are less sensitive +о moment orientation changes in the 
horizontal plane as compared to taa vertical plane. For 
examples, by keeping the vertical component of M constant at 
em angie of 50% down while varyirg the horizontal component 
Шгастіосі it was  fcund that ‘the horizontal component 
affected the positioning much less than variations іп the 
Enrucal. Illustrated in Figure 4.11, the x #s*imates of 
position improve as the horizontal component appro&ches the 
plazftcrm heading, while she y 2Stimates inprove as “һе hori- 
zontal component approaches alignment with magnetic North. 


Tho latter case is not shown. 


D. CHANGES IN SENSOR SPACING 


Changing the inter~sensor spacing had at most 2 weak 
EEG ct upon the position output, ven with svacing as short 
as 5 feet or as long as 200 feet. This wouid bs expected in 


a computer  simulaticn where sufficient precision is 


30 





шэг ос <рас хс 5 аза - аз 1220 Fast 5513 2:2%07 
became noticable, Тл ЕНЕ таит Та, senso 


Ћи «қаз «қыш» A 


о пцелпса *he localiza*ion estimates 


E. ALTITUDE CHANGES 


For a straight overhead pass, ров аса па п ғо жас оп 
appears best if a pass with а 500 зо 1000 2005 separation 
between the target and the E is made at CPA. At 
greater altitudes, СПЕ ЩЕ СМЕНА susters slightly with the 
increase in altitude, while at lesser ranges, usable posi- 
Каша па inferMation is obtained only reéelativsiy close to 
CPA. For exanple, for a vertical separation of 400 feet, 
positioning data was very inaccurazs at slant ranges greater 
than 1000 feet. 


Р. NON-OVERHEAD PASSES 


11 


This program gives poor results for passes that ars net 
mere yeove>: D2 directly oft to ons side cf she target. The 
HESS" ss*uqt2on observed in any variation is illustrated in 
pere 1.12 . To generate this figur 

that the actual x coordinate was 2pproximately equal 5с the 
actual v cocrdinate. The best pos H was achieved by 
overhead passes, but low off to the side passes car yield 


good results as well. 


6. AVERAGES 


The averages resulting from а number of runs are 
presented in Table I. These selected a 


e 
nor he variations discussed in {his Бов = пон 127714 


4 


iD 


С 
consistency. This is probably due in part to the existenc 
Ве Singularities. The variability is О 


1 
EBEN nts of zravel for ‘he runs. Nevertheless, for most 
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TABLE I 
Values of x and y Averaged over a Run 
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Figure 4.6 Platform Heading 1200 Magnetic. 
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Figure 4.7 Platform Heading 240° Magnetic. 
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Figure 4.9 Ш Vertically Up. 
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cases, hey d> show reasonable values, and Бү themseives іо 
Бао сен 11:55:31 Localization data. Ths problem for 110511945 
analysis would be to determine the bounds of a run in the 


absence of any other target positioning data. 


Н. IMPORTANCE OF CONSISTENCY 


An 2MpOrsant observation resulting  fzom ths analysis of 
the graphs, including those not shown as сеёхалрїез 11205115 


2 

Chapter, is tha- the rate of change or the calcu 

mates gives some indication of the valicizy of zhe posi- 
а 


Eson2rng data. If the x and y graphs have a sl 


сато, and Lf the zZ areph has a 51ора ОБ -459 for The calcu- 
lated coordinate, then positioning data is likely те be 
comparatively accurate. па Schzse атара отг: this 


Т 
condition siímul*aneously, then positioning date appears to 


Бе usable. 


1, IMPLEMENTATION 


Phe poss- Dle inplementa won oÍ the process described in 
Pees -п2525 would be to assume a magnetic moment for the 
Bert, then uz the sigral data through th= processirg 50 
generates positioning daca, and -hen compare these results to 
those cbtained by using the assum32i moment 32  gensza-3 a 
Signal via this simulation. if the agreement be-wesar results 
КОБЕ GOOG, chs could assume =hat one had assumed the correct 

= 


a 
gez moment; I$ not, one would assume a di 
n 


се ~nrouch the process again, until one had à reaso: 

terget moment, which would be ча 113 ас least until the 
target changed course. sanc a COMPU Ter -12== compaziscns 
coulé be made reasonably quickly, so this process mign- have 


apetat1ional possibilities. 
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Por the assumpticns made, Wemeosielouang estl2zma-es 
rsasonable. hiss, the averages taker of the three coor 
mates as calculated throughout the various runs show 
the x calculation averages are in generai at werst wit 
500 feet of the actual x value, whils the worst y calcu 
Босо а“ <=<2а025 nd to be within 1009 fetr of tie act 
value. The best calculated x and y values wars with: 
end 200 ==э= values. 


БЕ рес ЕТУ, O2 tie crue 
у 


ғ 
ends observed were as follows 


¡NTE s best to Ely che sensor EP m directly over 
the target, ООЛ О р е yelvy то Have а ‘distant’ 


ПЕ о аге 5o 5rla-fora з tion in the 
3 s The worst situation is wh 
cint that, az CPA, has a depressio 
р 


ect to the piatform. 


3. Bəs? results aze ортз іле =f headings ar? northerl 
ШІ soucheziy. Poor results are achieved with heading 
close to 609, 1200, 2400, and 3090. East-West headings 
eE noz as 1054 as North-South, do provide good posi 


МЕС па datas close to CPA. 


4. Good results та 


ћ а 
between 1000 and 2000 foot) CPA off +0 оле side while 
2 


а 
n the target is act 
Ч 


ү 


у Бә еєехрестзА4 lr there is a large 
Meerical monent, especially if it is vertically up. 
Poor positioning results are obtained with this simula- 
mom of the momen hes e large horizontal component. 


42 





Е в = оцена ре if the tazget monen*t is 
ШЕЕ р either the axis of flight or 
earth's field, although this effect is O 
by the more gross effects caused by che ct! 
A suggested apprcach for further stud 
noise into the Simulation in order to see h 
EN tore accuracy or the positicn calculations. E 
results remain reasonable, this approxi 
form of it may one day have real-wor 
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APPENDIX B 


COMPUTER PROGRAM FOR NUMERICAL OUTPUT 
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АРРЕМОТХ С 


EXAMELE OUTPUT OF APPENDIX B8 


Mecanico ehe ou put orka pogran in Appendix B is 


Shown below. The total output for a single cime 55:80 is 


SENSOR SLANT ЗЕМ ЕО CAL CULAT Ev 
NUMBER RANGE SIGNAL SIGNAL 
1 5650.0 0.00540-0% 0.5054с-0% 
2 ебе | 0e15560-C5 0„1612Е-03 
ES 5050.2 0.2682 0- Са 0.2 083Е- 05 
4 5099.9 д.5 9400-04 0.4 135Е- 04 
5 2(25:04:2 0,953660-04 0,9363С-09Ф 
ЕСКАТ ТОМ ТВОЕ APPROXIMA 
52-53 051321 -05 О ЗСО 
о ле со о О 0 Oeo! «8E-O 
ST zoe ГЕ 
52-51 0. 2312L-04 0,35314Е-0 
$3- са -,1258[-04 = ЕО 
в 9. 54261-94 02525525029 
EQNS X Y 
USED 2; POS РС5 
шта 3 ЏАС -1000,00 -4950.00 


123 ОЦ Рато ______=е0258 __=870483 24138141 


ИЕ Е == Five tows of data з0оүз ace ЁО 

"sensor" as ни mio, chant rsnge from 55325 sensor 20 
the target, The zeceilved signal as gensrated by equation 2.5 
and «he signal that weuld be generated usi С 


а 
Но Чесссірейі in СЦ ар:2т 3 “0 11 33:126 the prob 
NEX- six го 5 ә 


2 
between the recezlved and с 


соја = “(агат соправугов pur posas 
only) аа ои тз за under che htadiags 22 TRUS and 
APPROXIMATE tively. The £inal two lines sacw tha 


y respec 
Beeual ccerdinatss of the platforn relative t 


9 
соо аз атэ transfo-zamation hes not been dons az 


s 
Bema) as commared to the calculated position. These are the 
pers used by the program іп Appsndix A to plot the graphs 


Ex tne ссогазтате transformation is comple ec. 
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The seconi secticn of the output is shown 


beiow. 


EQUATION CUNS Х Ж 
МЈ МБЕ К | СОЕЕ Е СЕКЕТ COEFFICIENT USAR ICIENT 
1 УрФ =O. SUDt UZ -9.920%01 0.200ж21 
2 -Q0.2dutQ03 О. 13006 Os 130909 Ve EDU 4 
5 -0.бор+ -0.210%0і -0.,194%02 бе 3200: 
RATIO 
1/2 –0).132+01 -0., 22860253 -0.171:402 2% Эсет092 
1/3 -(428с%01 0. 44002 0.“9с%00 Da JOCUL 
2/3 9. 44Е+ 01 -0. 656701 -0.09с-92 UNE UL 
ШЕ 115: тігзе lines of the output represent the ccecfi- 
@sonts of the agua iors h. 6 = 3.83 in ths form of 
CONSTANT = X COEPF(x) + Y COEFF (y) + Z COEFF(z). In order чо 


determine the 


linear independence of these 


ratios of the coefficients are taces 
from the last three lines 
ОЕ 


cantly for 


be seen 
= па вопетелпс= апа 
SE 


аге 


ratios 
thrse equations, i- 


equations linearly 


throughout any given run, which 


(ша! time steps for that un. 


coefficients differ 
саа расой спаса 
independan 


nay bea 


equa-ion 
13 pairs. Sizce 15 па; 


оп ое 


holds 


ЕЕ ЕЯ byt pr ine ing 


Er This 


True 





APPENDIX D 
CRAMER'S RULE 


А. GENERAL 


Cramer's Rule, а method for solving systems of linear 
equations in two or more unknowns, sSstat2s that sach unknown 
can be expressed as the ratio of two determinates. The 


general form 9f a solution is shown below. 


Ax+By+Cz=0 


EX*Fy*t*Gz-H 


ix*JytKz-L 
BB C 

D = Eur 
т] 
CABE 


m ti 
t шщ 
N 
o 


B. EXAMPLE 


The follcwing example was taken from [Ref. 4]. Solve: 


Bia 2yt2Z=/ 





Let D be the denominator. 


1173 = 2: 12 | 
р = їз. 1, A = -17 
сты; 
i 2 20| 
| 6 1 1 | =з 
Хэ-г 022-1-2 | = -- 53 
----------- -17 
р 
ри 2 
|) 2-2 en 
IA ӘР 
D 
3 -2 7 | 
1 1 6 | - 17 
2 = {2-1 2 | а --- = 1 
----------- -17 
р 
CHSCK: 
3х-2у%22-7 
903) -2 (2) *2(1) =7 
9-4+2=7 
1=7 
КУ %2-6 
525146 
6-6 
2221-2252 
2(3)-2-2(1)=2 
2-2-2-2 
2-2 


С. SUBROUTINE "CRAMER'S ROLE" CHECK 


EMunserg the values in Section 8 through пе subroutines 
of the programs in Appendix А апі Aopendix B yielded the 


Sam? results as the examole. 
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